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£2 '. Abstract 

'NT . 

' Measuring the "^Be solar neutrino flux is crucial towards solving 

Q>^ . the solar neutrino puzzle. The Borexino experiment, and possibly the 

KaniLAND experiment, will be capable of studying the ^Be neutrinos 
Q^' in the near future. We discuss (1) how the seasonal variation of the 

• Borexino and KamLAND data can be used to measure the '''Be solar 

D , neutrino flux in a background independent way and (2) how anomalous 

seasonal variations might be used to discover vacuum neutrino oscil- 
lations, independent of the solar model and the measurement of the 
^ . background. In particular, we flnd that, after three years of Borexino 

^ [ or KamLAND running, vacuum neutrino oscillations can be either es- 

tablished or excluded for almost all values of (sin^ 26, Arn^) preferred 
by the Homestake, GALLEX, SAGE, and Super-Kamiokande data. 
We also discuss how well seasonal variations of the data can be used 
to measure (sin^ 29, Am"^) in the case of vacuum oscillations. 
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1 Introduction 



The question whether neutrinos have non-zero mass has been an outstand- 
ing issue in particle physics for many decades. Recently there have been 
new exciting developments in the indirect search for neutrino masses via 
neutrino oscillations. Major progress has been achieved by studying atmo- 
spheric neutrinos, culminating in the announcement of evidence for muon 
neutrino oscillations by the Super-Kamiokande collaboration The most 
striking signal presented in |I| is the up-down asymmetry of the atmospheric 
muon neutrino flux. The choice of this particular quantity eliminates many 
theoretical uncertainties and the final result is very robust. In fact, at present 
time, this result represents perhaps the best evidence for physics beyond the 
Standard Model. 

Another very active area of research is the study of neutrinos coming from 
the Sun. Ever since the Homestake experiment reported its first results, 
there has been disagreement between theoretical predictions and measure- 
ments of the solar neutrino flux. For many years, however, it was not possi- 
ble to determine if the observed discrepancy was due to problems with the 
experiment and/or with the modeling of the Sun, or if it was, in fact, a sign 
of new physics. In the last decade other neutrino experiments, Kamiokande 
0, GALLEX 0], SAGE 0, and more recently Super-Kamiokande 0, have 
also measured the solar neutrino flux, with different energy thresholds and 
using very different techniques. All four experiments confirm a deficit in the 
observed number of solar neutrino induced events. Moreover, it has recently 
become clear that it is virtually impossible to concoct a solar model which 
would fit all the data 0, |]. On the other hand, the results of all experiments 
can be explained by assuming that the electron neutrino oscillates into a 
different flavor state. 

There are two neutrino oscillation scenarios that are capable of faithfully 
explaining the solar neutrino data [0. One scenario makes use of the MSW 
effect [T^, where the electron neutrino conversion into another neutrino 
flavor is due to flavor dependent interactions with solar matter. The other 
scenario assumes that the neutrino oscillation length is comparable to the 
Earth-Sun distance, and simple vacuum oscillations are sufficient. This sce- 



nario is also known as the "just-so" solution [|TT|. Both scenarios allow the 
electron neutrinos to oscillate into other active species or sterile neutrinos. 
The solar neutrino energy spectrum is determined by several nuclear re- 



actions which take place in the Sun's core [O, and different experiments 
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are sensitive to neutrinos produced by different nuclear reactions. Super- 
Kamiokande, for example, a very large water Cherenkov detector, is currently 
sensitive to solar neutrinos with energies slightly above 5.5 MeV. Almost all 
neutrinos it detects come from the decay of (®B — >■ ^Be* + e"*" + z/e).Q 

Another solar reaction that gives rise to neutrinos is the process of elec- 
tron capture by a '^Be nucleus C^Be ^Li +z/e). Neutrinos from this 
reaction have energies below the Super-Kamiokande threshold, but are ac- 
cessible to the radiochemical experiments Homestake, GALLEX, and SAGE. 
If one naively assumes that the suppression in the neutrino flux is due to 
the suppression of individual neutrino sources (^B, ^Be, etc) in the Sun, the 
combination of the Super-Kamiokande data with that of the radiochemical 
experiments indicates that the flux of ^Be neutrinos is virtually absent M, IT^ 



(the best fit value of the ^Be flux is in fact negative!). In the case of neutrino 
oscillations, all solutions to the solar neutrino puzzle indicate that the ^Be 
neutrino flux is suppressed, in some cases very strongly. Thus, at present, 
there is great demand for experiments that would accurately measure the flux 
of the ^Be neutrinos. Two upcoming experiments, Borexino and KamLAND, 
may have the capability to do exactly that. 

In this paper, we present a quantitative study of what can be accom- 
plished by measuring the seasonal variations of the ^Be neutrino flux at 
Borexino and KamLAND. Seasonal variations of the solar neutrino flux are 
of course expected, due to the Earth's eccentric orbit. The number of neutri- 
nos of all flavors reaching the Earth is larger when the Earth is closer to the 
Sun than when it is farther away, and should vary as l/L^. In the case of no 
neutrino oscillations or of the MSW solution to the solar neutrino puzzle, the 
number of ^Be solar neutrino induced events is supposed to vary according 
to the law, following the variation of the total neutrino flux. This will 
be referred to as the "normal" seasonal variation. 

If vacuum oscillations are the solution to the solar neutrino puzzle, large, 
anomalous seasonal variations of the number of ''Be solar neutrino induced 



events might be detected 14]. It is well known that neutrino oscillation 



effects depend on the distance to the neutrino source, and different Earth- 



Sun distances may yield very different survival probabilities [15, 16 1. The 



anomalous seasonal variation effect should be more pronounced in ''Be neu- 
trinos than in ^B neutrinos (the latter was recently studied in fl^)- This is 



* There is a small fraction of the neutrinos that can be detected at Super-Kamiokande 
coming from the reaction '^He +p ^"^He +e+ + i/g. 
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due to one important feature which distinguishes ^Be neutrinos from and 
other abundant types of solar neutrinos: because they are produced as part of 
a two-body final state, the neutrino energy spectrum is mono-energetic.[| The 
details will become clear when we discuss the anomalous seasonal variation 
effect, in Sec. |[ 

In the case of no anomalous seasonal variations, if one has enough statis- 
tics and a small enough background, the time variation of the data can be 
used to measure the solar neutrino flux, given that the number of background 
events is constant in time.0 We will analyze how well Borexino and Kam- 
LAND can perform this type of measurement. We are particularly interested 
in analyzing the relevance of this technique when the number of electron 
neutrinos reaching the detector is very suppressed with respect to the Stan- 
dard Solar Model predictions, as might be the case if there are z/g — > z/^,r 
oscillations for the small angle MSW solution.^ 

The paper is organized as follows. In Sec. ^ we discuss how seasonal 
variations might be used to determine the solar neutrino flux at Borexino 
and KamLAND, in such a way that no separate measurement of the number 
of background events is required. In Sec. |^ we analyze the effect of the vacuum 
oscillation solutions to the solar neutrino puzzle on the annual variation of 
the number of detected events at Borexino and KamLAND. In particular 
we describe the region of the (sin^ 2^^, Am^) parameter space where vacuum 
oscillations can be discovered by studying the seasonal variations of the data. 
In Sees. ^ and ^ we describe how the measurement of the seasonal variation 
of the ^Be solar neutrino flux may be used to either measure the neutrino 
oscillation parameters, sin^ 26 and Am^, or exclude a large portion of the 
(sin^ 29, Am^) parameter space. In Sec. |] we summarize our results and 
conclude. 



^In fact there are two distinct neutrino energies, 0.383 and 0.862 MeV, corresponding 
to different final states of the ^Li nucleus. Borexino and KamLAND are only sensitive to 
the higher energy component. 

^Actually, a time-dependent background is also acceptable, as long as it can be moni- 
tored and understood well enough. 

§If Ve oscillates into sterile neutrinos, the suppression is even more pronounced, due to 
the absence of neutral current v^^r-e elastic scattering. We do not consider oscillations 
into sterile neutrinos in this paper. 
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2 Measuring the '^Be Solar Neutrino Flux 

As was already pointed out, measuring the flux of ''Be neutrinos is crucial 
towards understanding the solar neutrino puzzle. Borexino fl^ plans to do 
this measurement by using 300 tons of organic liquid scintillator to detect 
recoil electrons from elastic v-e scattering. Since the scintillator has no di- 
rectional information, and the signal is characterized only by the scintillation 
light produced by the recoil electron, the background has to be kept under 
control. This places a very stringent constraint on the radio-purity of the 
scintillator and on the activity of the material in the detector. Borexino 
anticipates 100 tons of fiducial volume for detecting solar neutrinos. 

KamLAND [|19], which was originally conceived as a reactor neutrino 



experiment with an unprecedented baseline (170 km on the average), may 
also be able to study ^Be solar neutrinos, if rigorous yet attainable require- 
ments on the radio-purity and activity are met. We assume throughout the 
paper that KamLAND will use 600 tons of fiducial volume for detecting so- 
lar neutrinos (the size of the fiducial volume will depend on the background 
rate, which is currently unknown). We concentrate our analysis on Borexino, 
which is an approved dedicated solar neutrino experiment, and discuss Kam- 
LAND, whose uses for solar neutrino studies are at present being proposed 



20|, as a possible higher statistics improvement. 

It is important to define what is meant by "measuring the ^Be solar 
neutrino flux." In reality, what the experiments are capable of measuring 
is the number of recoil electrons induced by solar neutrino interactions in a 
given recoil electron kinetic energy range (kinematic range). This information 
can only be converted into a solar neutrino flux measurement if one knows 
the flavor composition of the solar neutrinos Explicitly, assuming that 
the solar neutrino flux is composed of z/g (with fraction P) and z/^^r (with 
fraction Q = 1 — P) , 

#recoil elect rons/time = $ x {Pa^^^-e + (1 - P)ai,^^-e)A^e, (2.1) 

where $ is the neutrino flux, N^. is the number of target electrons, and 



a.,-e = / dT[ — \ , (2.2) 




with ( ^ j being the differential cross section for v^-g scattering for a given 
kinetic energy T of the recoil electron. Tmin and Tmax define the kinematic 
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range. In the case of neutrino oscillations, P is the the survival probability 
for electron neutrinos, while 1 — P is the probability that will oscillate into 

If the flavor composition of the flux is not known, all that can be quoted 
is the effective neutrino flux, $eff, which is calculated from the number of 
measured recoil electrons assuming that there are only electron neutrinos 
coming from the Sun. Explicitly, 

^ #recoil electrons/time ^ ^ /p ^ _ -_^\ ^3 3) 

Clearly, if P = 1, $eff = It is important to remember that cXi/^^-e/ci/e-e < 1 
and therefore $eff < ^- The ratio of the neutrino elastic cross sections de- 
pends on the energy of the incoming neutrino and the kinematic range to 
which each particular experiment is sensitive. For = 0.862 MeV and 
the Borexino (KamLAND) kinematic range 250-800 keV (280-800 keV), 
'^vij.,T-e/ <^ve-e = 0.213(0.214). It is this effective electron neutrino flux, $eff) 
that is referred to, throughout this paper (and in general), as the ^Be solar 
neutrino flux. 

In order to determine the number of recoil electrons induced by solar 
neutrino interactions, it is crucial to determine the number of background 
events. The number of background events can be estimated by various tech- 
niques, which we do not address in this paper. It is worthwhile to point out, 
however, that this is a very difficult process and it would be highly desirable 
to have an independent way to determine the ^Be solar neutrino flux in order 
to make the final results more convincing. This may be possible if one looks 
at the seasonal variation of the number of detected events. 

In the following, we study the seasonal variation of the event rate as a 
means to measure the ^Be solar neutrino flux. The distance between the 
Earth and the Sun varies slightly over seasons because of the eccentricity 
of the Earth's orbit. The perihelion (when the Earth is closest to the Sun) 
occurs around January first. The eccentricity of the Earth's orbit is e = 0.017, 
and hence the distance varies as 

L = Lo(l - ecos(27ri/year)), (2.4) 

to the first order in e. Here, t is the time measured in years from the per- 
ihelion, and Lq = 1.496 x 10^ km is one astronomical unit. The neutrino 
flux varies as 1/L^ and hence shows a seasonal variation of about 7% from 
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minimum to maximum. The change in the Earth-Sun distance between the 
aphehon and the perihehon is given by 



AL = Linax - ^min = 2eLo = 5.1 X 10^ km. (2.5) 

By fitting the event rate to the seasonal variation expected due to the eccen- 
tricity, 

B + S(^)\ (2.6) 



L 

one can extract the background event rate B and the signal event rate 5* 
independently. As long as the detector is monitored well and its performance 
is sufficiently stable, this method will be only limited by statistics. 

Borexino expects 53 events/day^] according to the BP95 |2^ Standard 



Solar Model (SSM), together with 19 background events/day ||rH|, after the 
statistical subtraction of the known background sources. This is done by 
pulse shape discrimination against the a-particle background and the mea- 
surement of Bi-Po pairs via a-jS coincidence. This in turn allows the sta- 
tistical subtraction of processes in the ^ss-y g^^^^ 232rpj^ chains which are in 
equilibrium. It is also assumed that the experiment can achieve a radio- 
purity of 10-^6g/g for U/Th, 10~i^g/g for ^°K, ^^C/^^C = 10-^^ and no Rn 
diffusion. For KamLAND we use 466 events/kt/day for the signal and 217 
events/kt/day pO| for the background under similar assumptions but with 
larger cosmogenic background (especially ^^C) and some Rn diffusion. As- 
suming 600 t of fiducial volume, we expect 280 signal events/day and 130 
background events/day. Throughout the paper, we will assume that the num- 
ber of background events is either constant in time or its time dependence 
is sufficiently well understood by monitoring. We neglect systematic effects 
and assume that there are only statistical uncertainties. 

Under these assumptions. Fig. |l| depicts a simulation of the seasonal vari- 
ation of the "data" for both Borexino and KamLAND, after three years of 
running. The plots are for the case of the small angle MSW solution to the 
solar neutrino puzzle, where the z/g's produced by ^Be electron capture inside 
the Sun have almost completely oscillated into or i>r, and the event rate is 
reduced to 21.3% (21.4%) of the SSM prediction at Borexino (KamLAND). 

^For simplicity, we neglect the contribution of solar neutrino sources other than ^Be 
electron capture throughout the paper. In particular we neglect the contribution of neu- 
trinos produced in the CNO cycle, which is about 10% of that from the ^Be neutrinos. 
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Figure 1: The simulated seasonal variation of the ^Be flux for the case of the 
small angle MSW solution, for three years of Borexino (left) and KamLAND 
(right) running. The inset shows the measured flux of ^Be neutrinos from 
the flt to the seasonal variation of the event rate (point with error bar) and 
the SSM prediction (shaded band). 



In the flt to the data, both the background and the Be flux are allowed to 
float. 

This analysis can be repeated for different values of the ^Be flux, or, 
equivalently, for different survival probabilities for Uf.- Fig. |^ depicts the 
expected 1 a statistical accuracy of the ''Be flux measurement, together with 
the central value normalized by the SSM prediction, as a function of the 
survival probability for z/g. We emphasize that this measurement technique 
assumes no knowledge of the background. 

The important information one should obtain from this analysis is if one 
can indeed measure a nonzero ^Be solar neutrino flux. For example, in the 
case of the small angle MSW solution, the z/g survival probability is very close 
to zero and, assuming the expected number of background events, Borexino's 
measured neutrino flux is less than 1.5 a away from zero. The situation at 
KamLAND is much better, and in the case of the small angle MSW solution 
a healthy 3 sigma-away-from-zero measurement of the flux is obtained, if 
the background is as low as expected. The signiflcance of the measured flux 
increases for larger survival probabilities, as in the case of the large angle 
and the low Am? MSW solutions. 

A similar analysis can be performed in order to determine how many 
background events each experiment can tolerate in order to claim a solar 
neutrino flux measurement which is 3 a away from zero. Fig. ^ depicts the 
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Figure 2: The expected 1 a statistical accuracy of the ^Be neutrino flux mea- 
surement, together with the central value normalized by the flux predicted 
by the SSM, as a function of the v^. survival probability at Borexino (left) 
and KamLAND (right), after three years of data taking. 
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Figure 3: The maximum number of background events allowed per day at 
Borexino (left) or KamLAND (right), for 3 years of running, in order to 
measure a solar neutrino flux which is 3 cr away from zero. The dashed lines 
indicate the currently anticipated number of background events per day. 
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maximum number of background events per day allowed for 3 years of Borex- 
ino or KamLAND running. It is worthwhile to comment that, in the case 
of Borexino and the small angle MSW solution (P ~ 0), a 3 sigma-away- 
from-zero measurement of the neutrino flux is not attainable in three years, 
even in the case of no background (note that for P ^ 0.05 the required max- 
imum background to achieve a three a measurement of the flux is negative, 
i.e., impossible to achieve). Therefore, for Borexino, this simple, background 
independent analysis using the seasonal variation of the data is not particu- 
larly powerful in the case of the small angle MSW solution, due to statistical 
limitations. 

3 Sensitivity to Vacuum Oscillations 

In this section we study the discovery potential of the Borexino and Kam- 
LAND experiments in the region of Am^ corresponding to the vacuum os- 
cillation solution to the solar neutrino problem. In this case, the pattern of 
seasonal variations can be very distinct from the normal pattern discussed 
in the previous section. 

The basic idea is the following. The survival probability P for an electron 
neutrino in the case of neutrino vacuum oscillations between two flavor states[| 
is given by 

P = 1 - sin2 2^sin2 (^1.27Am2^^ , (3.1) 

where the neutrino energy E is in GeV, the distance L in km, and the differ- 
ence of masses- squared Am^ in eV^. Model-independent analyses of all solar 
neutrino data show the need for an energy-dependent suppression of the 
flux. The "just-so" solution achieves this by choosing Am? such that the 
corresponding neutrino oscillation length 

i^osc = 7 — ^ = 2.47 X 10^ km x — - — - — 3.2 

1.27Am2 VlOMeV/V / 

is of the order of one Astronomical Unit (1 a.u. = 1.496 x 10^ km); hence 
the name "just-so" . More specifically, the oscillation length is assumed com- 
parable to 1 a.u. for ^B neutrinos [Ey ~ 10 MeV); at the same time, the 

*One can assume the more complicated case of oscillations between three neutrino 
flavor states. In this paper we limit our studies to the case of oscillations between two 
flavor states. 
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Figure 4: Illustration of the effect of vacuum oscillations on the shape of the 
seasonal variation of the solar neutrino data. The points with statistical error 
bars represent the number of events/month expected at Borexino after 3 years 
of running for /S.rn? = 3 x 10~^° eV^, sin^ 26 = 1. The histogram in (a) shows 
the number of events predicted by the SSM without neutrino oscillations, plus 
the number of anticipated background events. The histogram in (b) shows 
the same quantity after adjusting the solar neutrino flux and the background 
rate so as to minimize the value of x^, as explained in the text. The difference 
between the case with oscillations and the one without oscillations is still 
apparent. 



oscillation length of ''Be neutrinos {E^, = 0.862 MeV) is an order of mag- 
nitude smaller and, for sufficiently large Am^, can be comparable to the 
seasonal variation of the Earth-Sun distance due to the eccentricity of the 
Earth's orbit, AL (see Eq. ( p.5|) ). As a consequence, the flux of '^Be neutrinos 
detected on the Earth may exhibit an anomalous seasonal variation, beyond 
the normal effect discussed in the previous section. 

Such anomalous variation could serve as a unique signature of vacuum 



oscillations [0, |14| . Moreover, as we will show in this section, both Borexino 
and KamLAND will be able to cover a large portion of the "just-so" param- 
eter space, even without relying on a particular solar model or estimate of 
the background rate, just by analyzing the shape of their data. In this sense 
the discovery of an anomalous seasonal variation at one of these experiments 
would be as robust a result as the Super-Kamiokande measurement of the 
up-down asymmetry for the atmospheric muon neutrinos. 

To illustrate the main idea, we choose a particular point (Am^ = 3 x 
10^^° eV^, sin^ 2^ = 1) in the allowed region of the "just-so" parameter 
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spaceQ and compute the corresponding seasonal distribution of the neutrino 
events at Borexino after 3 years of running. We use the number of background 
events and the expected number of signal events (before the effect of neutrino 
oscillations) quoted in Sec. ^. The results are shown in Fig. ^ by the set of 
"data" points with error bars; each point represents the number of events 
expected in a given month and the vertical error bars show the corresponding 
statistical uncertainties. The histogram in Fig. ^(a) shows "theoretical" event 
rates expected for non-oscillating neutrinos, provided the background rate is 
known accurately and the SSM prediction for the neutrino flux is trusted. 
One can see that under these assumptions vacuum neutrino oscillations with 
Am^ = 3 X 10^^° eV^, sin^ 29 = 1 would be trivial to discover. 

More importantly, the experiment would be able to claim the discovery 
even without relying on an estimate of the background rate or the value of 
incoming neutrino flux predicted by the SSM. It is intuitively obvious from 
the figure that the vacuum oscillation "data" points cannot be fit by the 
"theoretical" curve even if the background and the solar neutrino flux are 
varied freely, unless one assumes neutrino oscillations. This can be quantified 
as follows. For a given background rate b and signal event rate s, we define 
the value of the fit for an "average" experiment: 

x\s, b) = iVd.o., + j: (d^-b-s.h,r ^ 

where A^bins is the number of bins, A^d.o.f. is the number of degrees of freedom, 
di is the average expected number of neutrino events in the ith bin, and hi 
is given by hi = //_i(l — e cos(27rx/A^bins))^'^a; . The constant term A^d.o.f. in 
Eq. (|37^ ) is added to take into account the effect of statistical fluctuations in 
the data. In a single experiment, statistical fluctuations make the number of 
neutrino events in the ith bin slightly different from di, and is computed 
by an expression similar to Eq. (|3.3|), with di replaced by the number of 
events measured in the ith bin and without the constant term, A^d.o.f.- In 
our analysis, however, we are interested in the sensitivity of an "average" 
experiment. As proven in Appendix |^, averaging over many experiments 
results in the definition of given in Eq. ( ^.3[ ), with the constant term Ad.o.f.- 
This agrees with the conventional wisdom that, if a function describes data 
correctly, the average expected value of should be equal to the number of 

^Based on the analysis of the total rates in the Homestake, GALLEX, SAGE, and 
Super-Kamiokande experiments. See Fig. 5 in H. 
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degrees of freedom. Given this definition, we can choose values of s and b 
that minimize the x^] the only restriction imposed is that both s and b be 
non-negative. For the case at hand the minimum occurs when b is zero and s 
is 0.95 times the SSM prediction (see Fig. |(b)). As expected, even after this 
change the "data" points and the histogram are very different. (Numerically, 
= 2935 which for 10 degrees of freedom implies a confidence level of 
1 - 9 X 10-626!). g 

We now extend this approach, and scan the entire (sin^ 29, Am^) plane 
(for an earlier work with a more simplified analysis which does not consider 
the presence of background, see |15[)- In the analysis below, we follow the 
same steps as before: the "data" is simulated according to the expected num- 
ber of background and signal events, plus the effect of neutrino oscillations, 
for each value of (sin^ 29, Am^), binned into a certain number of bins Ni^ms, 
and then compared to the "theoretical" predictions for the case of no os- 
cillations. The is computed according to Eq. ( p. 3D and minimized with 
respect to both the signal (s) and background (6). The confidence level (CL) 
corresponding to the minimal value of and iVd.o.f. = A^bins — 2 degrees of 
freedom is then determined, and the region in which the CL is less than a 
given number is isolated. This case, when both the number of background 
events and the incoming solar neutrino flux are considered unknown in the 
fit, is the most conservative one, and yields the smallest sensitivity region. 
Later we also study less conservative cases, where we assume in the "data" 
analysis that the incoming neutrino flux is the one predicted by the SSM 
and/or that the background rate is known. 

We now apply this most conservative procedure to study the experimental 
reach of Borexino after 3 years of operation. In Figure ^ we show the results 
of the scan for 95% and 5 a CL. As one can see from the figure, even at 
5 a CL a large portion of the parameter space above Am^ ~ 10~^° eV^ is 
covered (white region). In this region the neutrino oscillation length Lose is 
smaller than the seasonal variation of the Earth-Sun distance AL. On the 
other hand, below Am^ ~ 10^^° eV^ one can see a series of spikes protruding 
through the sensitivity region. It is important to understand the origin of 
these spikes. Since we adjust the level of signal and background in the fit, we 
are not sensitive to the absolute event rate, only to its variation during the 
year. For Am^ < 10~^° eV^ the oscillation length is larger than AL and the 

■l-This number is, of course, unrealistic, and the true confidence level in this case will 
be dominated by systematic effects. 
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Figure 5: The sensitivity region of the Borexino experiment in 3 years, if 
the analysis does not assume any knowledge of the background rate or the 
incoming solar neutrino flux. In the unshaded region the "data" is at least 
5 a away from the best no-oscillations fit. In the lightly shaded region the 
discrepancy is greater than 95% CL but less than 5 a CL. 

amplitude of the variation of the event rate is roughly proportional to the first 
derivative of Eq. (|3.1| ) with respect to L. In the regions where this derivative 
nearly vanishes, the amplitude of the variations is small and the signal is 
indistinguishable from the case of no oscillations. This explains why the loss 
of sensitivity occurs not only when the neutrinos undergo approximately an 
integer number of oscillations as they travel to the Earth (Am^ = nx 0.143 x 
10~^° eV^), but also when the number of oscillations is close to a half-integer 
(Am^ = (n + 1/2) X 0.143 x 10"^° eV^). In the latter case the absolute 
neutrino flux is maximally suppressed, but the magnitude of the seasonal 
variation is small. | 

Given this explanation, one would expect that the spikes corresponding to 

§ Notice that the regions preferred from the global fits have the absolute ^Be neutrino 
flux suppressed. See Figs. || and |l^. 
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a half-integer number of oscillations should become shorter if in the analysis 
we choose to rely on the SSM prediction of the incoming neutrino flux and/or 
on the anticipated background rate. It is straightforward to incorporate the 
knowledge of both quantities and their uncertainties in our procedure. For 
example, to impose the value of the incoming neutrino flux predicted by the 
SSM, we modify the expression of in E<1- (p.3|) by adding an extra term: 



x'{s,b)^X%s,b)+^-^^, (3.4) 

so 



where s and b are the values of the signal and background with respect to 
which we later minimize x^? ■So is the SSM prediction for the signal, and 
ago is the uncertainty in sq. The rest of the analysis is carried out un- 
changed, except that the number of degrees of freedom is increased by one 
to Nd.o.f. = Nuns — 1- To use both the incoming flux predicted by the SSM 
and the anticipated background rate, two terms are added to Eq. ( p.3|) and 
the number of degrees of freedom is increased by two to iVd.o.f. = iVbins- 

The results of the calculation are shown in Fig. |^. The uncertainty on 
the solar model prediction of the ''Be neutrino flux is taken to be 9% P3 



while the uncertainty on the background is 10% pO[. As expected, the odd- 



numbered spikes do become shorter. The one possibility not shown in the 
plot is the situation when one only assumes knowledge of the background 
rate. In this case the spikes become significantly thinner, although their 
length remains virtually unchanged. 

In order to extend this analysis to values of Am^ > 10~^ eV^, several 
issues must be confronted. We will next address these issues one by one, and 
illustrate the discussion in Fig. ^ 

The first and the most obvious point is that the number of bins needs 
to be changed. The reason is that the frequency of the seasonal variations 
increases with Am^, and above some value (Am^ ~ 8 x 10^^*^ eV^, for 12 
bins) integration over the bin size washes out the effect. To avoid this, we 
change the number of bins from 12 to 365. After the change, the effect of 
binning kicks in at Am^ ~ 2.4 x 10~^ eV^, as curve 1 in Fig. |^ illustrates. 

Next, there are two physical effects one must take into account: one is 
the interaction of the neutrinos with solar matter (the MSW effect), and the 
other is the finite width of the ^Be solar neutrino line. One may worry about 
the wash-out of the seasonal variation effect due to the finite size of Sun's 
core. However, matter effects make the core size effect irrelevant because the 
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Figure 6: The sensitivity reach of the Borexino experiment after 3 years 
of running (at 95% confidence level). The three cases considered are: no 
knowledge of either the background rate or the incoming solar neutrino flux 
(the covered region is white); assumption that the incoming solar neutrino 
flux is the one predicted by the SSM, with 9% uncertainty (the covered region 
is white + light gray); assumption that the background rate is known with 
10% uncertainty and the incoming neutrino flux agrees with the SSM, with 
9% uncertainty (the covered region is white + light gray + medium gray). 



mixing angle in the Sun's core is small and the oscillations effectively start 
at the level-crossing point (see Eq. (|3.6| )).P1 

When a z/g is created by the electron capture process in the core of the Sun, 
its Hamiltonian is dominated by the matter effect ^/2Gpnf, [n^, is the electron 
number density) if Am^ -C 10~^ eV^ for ^Be neutrinos. We restrict ourselves 
to Am^ < 10"'' eV^ in the following discussions, as the final sensitivity due 
to the anomalous seasonal variation is limited by ^ 10~® eV^ as will be 
seen later in this section. Then the mass mixing effect can be completely 



^We thank E. Lisi and L. Wolfenstein for pointing this out to us. For earher papers on 



this particular point, see ||9|], |11 , and in particular, 
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Figure 7: The relative roles of the binning effect, the linewidth effect, and 
the matter effect, as explained in the text. 

ignored at the time of the neutrino production, and one can safely take 
the produced neutrino to be in a Hamiltonian eigenstate (the one which 
corresponds to the larger energy in the Sun's core). As it propagates through 
the Sun, the neutrino follows the instantaneous Hamiltonian eigenstate (in 
the adiabatic approximation), and exits in the heavier mass eigenstate, z/2 = 
UgSinO + UnCosO. It also has a finite amplitude A^. for hopping to the other 
Hamiltonian eigenstate. The neutrino state that exits the Sun can therefore 
be written as 

i^exit = ^cJ^i + Bciy2, (3.5) 
with the unitarity constraint l^d^ + |i?cP = 1- Out of the Sun, the two mass 
eigenstates develop different phases due to the mass difference, e-«^'"^*/2E^ _ 
Therefore the neutrino state that arrives at the Earth is given by 

i^axrival = A,U^ + S.Z.^e"^^'"'^/'^'' , (3.6) 

up to an overall phase factor. The distance L is between the point of level 
crossing and the Earth. Finally, the survival probability of the electron 
neutrino is determined by the z/g component of i^arrivai, and hence 
P = |A,cos^ + 5eSin^e-'^'"'^/2E.|2 

= \A^\^cos^e + \B^\^sm^e + 2ReAlB^e-'^'^'^/^^-' sine cose. (3.7) 
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Figure 8: The contour plot of the hopping probabihty Pc = 0.1, 0.2, . . ., 0.9, 
for the ''Be neutrino energy, using the exponential-profile approximation for 
the electron number density and Eq. ( ^.9| ). 



Since l-BcP is the hopping probability between two Hamiltonian eigenstates 
in the Sun Pc, one can rewrite the formula using Pc and an additional phase 
factor A*Bc = J Pdl - Pc)e-'\ 



P = Pc cos^ e+{l-Pc) sin^ e + 2JPc{l-Pc) sin 6 cos 6 cos I^HLLl + ^ | . 

(3.8) 



An approximate formula for Pc was given in [25| using the exponential density 
profile of the Sun, 

(3.9) 



1 - e-^ 
with 

, = 2.r„^ = 1.22 , (3,10) 

where we consider the exponential-profile approximation for the electron 
number density in the Sun Ue (x exp(— r/ro), with tq = i?0/lO.54 = 6.60 x 
10^ km, given in ||26|. Fig. |^ shows the contours of Pc on the (sin^2^, Am^) 
plane for the ^Be neutrino energy E^, = 0.862 MeV. 
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The most important consequence of the matter effect is that the vacuum 
oscillation is suppressed when ~^ (adiabatic limit). The origin of the 
suppression is simple. When Pc is small, the neutrino state that exits the Sun 
is nearly a pure V2 state. Since it is a mass eigenstate, only its phase evolves 
in time and no oscillations take place. The Ve survival probability then is 
simply given by the content of z/2, which is nothing but sin^ 6', without 
anomalous seasonal variations. Therefore, the sensitivity to the anomalous 
seasonal variation is reduced in the region with small Pc- When Am^ is 
small, on the other hand, the situation is in the extreme non-adiabatic limit, 
and Pc — > cos^ 6^. Then Eq. ( p.8| ) reduces to Eq. As Am^ increases, Pc 

becomes smaller than cos^ 9, which enhances the vacuum oscillation effect in 
the small mixing angle region. Curve 2 in Fig. |^ includes the matter effect 
and indeed indicates a reduced sensitivity for large sin^ 26 (small Pc) and an 
enhanced sensitivity for small sin^ 26 (where Pc starts deviating from cos^ 6)^ 

The second effect is the finite width of the ^Be line. To give some prelim- 
inary idea about the relative size of this effect, we first consider a simplified 
model. We assume for a moment that the only source of the line broaden- 
ing is the Doppler shift of neutrino energies arising from the thermal motion 
of the ^Be nuclei. Since the energy is shifted to — > £^(1 -f- Vz/c) and 
the probability distribution of the velocity along the line of sight Vz is pro- 
portional to exp(— mf^/2/cT), the resulting line profile will be a Gaussian 
exp(— mc^(£' — Eq)"^ / {2kT El)) . Taking the temperature to be 15.6 million 
Kelvin (the temperature in the center of the Sun) and integrating over the 
line profile, we obtain curve 3 in Fig. |^. The sensitivity loss now occurs at 
Am? ~ 1 X 10~® eV^, demonstrating that this effect is more important than 
the matter effect. 

This naive model is actually incomplete; there exists another very impor- 
tant source of line broadening. Because the incoming electron in the process 
''Be + e~ — > ^Li -f- i/g has nonzero thermal kinetic energy, the center of mass 
energy of the reaction is greater than the one measured in the laboratory, and 
so the neutrino has a greater energy. The phase space distribution of elec- 
trons is governed by the Maxwellian factor ex]i{—Ec-/kT). This distribution 
has to be multiplied by the energy-dependent cross section, integrated over 
the phase space, and finally convoluted with the Gaussian arising from the 
Doppler effect. The resulting line shape becomes asymmetric, with a Gaus- 

II In the numerical scan, we ignored the additional phase factor S, because its effects are 
neghgible p^ . 
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sian profile on the left (due to the Doppler effect) and an exponential tail on 
the right (due to the Maxwellian distribution of the electron energy). The 



issue was studied in detail in |^^, where the precise form of the profile was 
computed.|^ Repeating the calculation with this profile we generate curve 4 
in Fig. 1^. 

One can see that for this curve the cut-off occurs at smaller Am^. This 
behavior is expected, because the linewidth is now greater than when only 
the Doppler effect was included (curve 3 in Fig. 0). It is also worth noting 
that the cut-off sets in more gradually. This feature can be understood 
analytically by considering the Fourier transform of an exponential tail vs. a 
Gaussian tail. The details can be found in Appendix p. 

Finally, we can combine both the linewidth and the matter effects. The 
result is curve 5 in Fig. ^ As expected, the inclusion of the matter effect on 
top of the linewidth effect introduces only a small distortion to the sensitivity 
region. It is important to note that for Am^ ^ 5 x 10~^° eV^ none of the 
physical effects mentioned above affect the sensitivity region (curve 1 versus 
curve 5, in Fig. 0). 

We need to consider one last ingredient in the analysis. We again return 
to the issue of the number of bins. While choosing more bins is necessary 
for larger values of Am^, it simultaneously leads to a loss of sensitivity for 
smaller Am?. A better procedure is to use an optimum number of bins N^pt 
for each Am?. It can be shown that for our method of analysis (minimizing 

by varying the signal and background) and sufficiently large Am^ an 
approximate formula holds: Nopt — 2 x 10^°(Am^/l eV^). Of course, this 
formula should not be used when the optimal number of bins it predicts is 
too small. We choose to use 12 bins for Am^ < 6 x 10^^° eV^ and a variable 
number of bins A^'bins = 2 x W^{Am'^/l eV^) for Am^ > 6 x 10~i° eV^.^ 

In Fig. 1^ we show the entire sensitivity reach of Borexino after three years 
of running. The unshaded region will be covered at least at 95% CL, if in the 
analysis one allows the background and the incoming solar neutrino flux to 
float. The dark shading marks the additional portion of the parameter space 

**It turns out that other effects, such as collisional line broadening ||2^ or gravitational 
energy shift pTf , are unimportant. 

^^An alternative technique, which can be considered more rigorous but which would also 
be more computer intensive, is to Fourier transform the simulated data for every value 
of (sin^ 20, /S.w?) in the scan. One can then compare the intensities of the harmonics to 
those expected for the case of no oscillations. A description of this method can be found 
in |2g|. For our purposes varying the number of bins is sufficient. 
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that will be covered at least at 95% CL, if in the analysis one assumes both 
the anticipated background rate (10% uncertainty) and the SSM prediction 
of the ^Be solar neutrino flux (9% uncertainty). For Am^ ^ 5 x 10~^ eV^, 
the sensitivity to the anomalous seasonal variation gets lost because of the 
smearing due to the linewidth effect. However, there is an overall suppression 
of the flux due to the MSW effect in this region. To be sensitive to this overall 
suppression, we should return to a smaller number of bins to enhance the 
statistical accuracy. We therefore use 12 bins in this region.^ 

For comparison, we also superimpose the "just-so" preferred regions ob- 
tained by analyzing the total event rates in the Homestake, GALLEX, SAGE, 
and Super-Kamiokande experiments (Fig. 5 in [^). The plot shows that 
Borexino will be sensitive to almost all of the preferred region, even with- 
out relying on the SSM prediction of the incoming neutrino flux or on the 
knowledge of the background rate. Only two thin spikes protrude through 
the lower "islands". This overlap disappears completely when the antici- 
pated background rate and the SSM prediction for the incoming neutrino 
flux are used in the "data" analysis, in which case the entire preferred region 
is covered. 

Fig. |10| contains a similar plot for three years of KamLAND running. Be- 
cause KamLAND will have more statistics, it will be sensitive at 95% CL to 
the entire preferred region without relying in the analysis on the SSM pre- 
diction of the incoming neutrino flux or on the knowledge of the background 
rate. 

As mentioned earlier, the sensitivity to anomalous seasonal variations is 
completely lost for Am^ ^ 10~^ eV^. In this case the seasonal variation of 
the data is consistent with an average suppression of the incoming neutrino 
flux. In particular, in the case of the MSW solutions (10"^ eV^ < Am^ < 
10~^ eV^), no anomalous seasonal variations can be detected, as was implic- 
itly assumed in Sec. |[ 

At last, it is worth mentioning that the experiments will still be sensitive 
to a significant part of the preferred region even if the background rate or 
the incoming ''Be neutrino flux (for all flavors) turns out to be significantly 
different. For example, if the background rate at Borexino (KamLAND) 
turns out to be 30 (100) times higher than expected, the part of the preferred 

MOne can cover a slightly larger portion of the parameter space by using yet fewer bins. 
We chose 12 bins such that one can still verify the expected behavior of the signal 

even with a reduced flux, as we discussed in Sec. |^. 
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Figure 9: The final sensitivity plot for three years of Borexino running, after 
the inclusion of all effects limiting the reach of the experiment for large Am^. 
The white region corresponds to the sensitivity at more than 95% confidence 
level with both the incoming neutrino flux and background rate assumed to 
be unknown, and the dark region to the additional coverage when the SSM 
''Be flux and the background rate estimated elsewhere are used. Also shown 
are the regions preferred by the analysis of the total rates in the Homestake, 
GALLEX, SAGE, and Super-Kamiokande experiments 0. 
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Figure 10: The same as Fig|^, but for three years of KamLAND running. 
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region with Am^ > 10~^° eV^ will still be within the reach of the experiment, 
after three years of running. The sensitivity will be completely lost only if 
the background rate turns out to be three (four) orders of magnitude higher 
than anticipated at Borexino (KamLAND). The consequences of a ^Be solar 
neutrino flux smaller than predicted by the SSM can also be studied. If the 
^Be neutrino flux is for some reason suppressed by a factor of 5, KamLAND 
is still sensitive to the part of the preferred region with Am^ > 10~^° eV^, 
after 3 years of running. 



4 Measuring the Oscillation Parameters 

In this section, we address the issue of how well the two-neutrino oscillation 
parameters, sin^ 26 and Am^, can be extracted if the data collected at future 
solar neutrino experiments exhibits an anomalous seasonal variation. In order 
to do this, we simulate "data", according to the procedure developed in 
Sec. ^, for two distinct points in the parameter space, sin^ 20 = 0.7, Am^ = 
8 X 10"" eV^ ("low point") and sin^ 26 = 0.9, Am^ = 4.5 x 10-^° eV^ ("high 
point"). The low point is close to the best fit point presented in [0, while the 
high point is close to the point preferred by the Super-Kamiokande analysis 
of the recoil electron energy spectrum [RO]. The data is binned into months 



(12 bins per years), and Fig. |TT| depicts the annual variations for both the 
high and the low points, assuming three years of Borexino running. The 
no-oscillation case is also shown. 

In order to measure the oscillation parameters, we perform a 4 parameter 
(s, b, sin^ 29, and Am^) fit to the "data". The fit is performed by minimizing 

with respect to the incoming neutrino flux (s) and the background rate 
(6), as in Sec. ^, and computing it for fixed sin^ 29 and Am^. Fig. depicts 
the values of (sin^ 29, Am^) and the 95% CL contours (for two degrees of 
freedom), extracted from the "data" consistent with the low (light) and high 
(dark) points. Note that this is very different from what was done in the 
previous section. There, for each point in the (sin^ 26*, Am^) plane there was a 
different "data" set, and the "data" was fitted by a non-oscillation theoretical 
function. Here the "data" is fixed (either the low or the high point), and is 
fitted by a theoretical function which assumes neutrino oscillations. 

One should easily note that the extracted 95% CL contour for the high 
point consists of only two "islands", while for the low point one extracts a 
collection of "islands". The reason for this is simple. When Am^ ~ few x 
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Figure 11: Number of recoil electrons detected in a given month, for the low 
point, the high point (see text for description) and the case of no neutrino 
oscillations, after three years of Borexino running. 

10^^° eV^, the oscillation length is slightly smaller than AL (see Eq. ( p.5|) ). 
This means that the seasonal variation of the "data" has a very particular 



shape (as one may easily confirm by looking at Figs. ^, |Tl]), which cannot 
be easily mimicked by other values of Am^, even when the background rate 
and the incoming flux are varied in the fit procedure. 

When Am^ ~ several x 10~^^ eV^, the oscillation length is larger than 
AL, and the effect of seasonal variations is less pronounced. There is a 
collection of Am^'s that yields the same qualitative behavior. Because our 
fit procedure allows for the background rate and the neutrino flux to float 
freely, a good agreement with the "data" is met for a large portion of the 
parameter space. In order to make this discussion clearer, it is useful to 
describe in detail what happens to the number of electron neutrinos reaching 
the detector as a function of time. 

In the case of the low point: initially, when the Earth is at the perihelion, 
the Vf. survival probability is small and, as time progresses, monotonically 
increases until the Earth reaches the aphelion (after six months). The pro- 
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Figure 12: Measurement of the neutrino oscillations parameters sin^ 29 and 
Am^, assuming no knowledge of the SSM and the number of background 
events. The regions represent the 95% confidence level contours, for data 
consistent with the high (dark) and low points (light). The input points are 
indicated in the figure by the two crosses. See text for details. We assume 3 
years of Borexino running. 

cess happens in reverse order in the next six months, as expected. There 
are many other values of the oscillation length, i.e. Am^, such that the sur- 
vival probability monotonically increases for increasing Earth-Sun distance 
and therefore a similar qualitative behavior is to be expected. The main 
quantitative difference is in the ratio of the number of events detected in the 
perihelion and in the aphelion, which may be accounted for by varying the 
background rate and the incoming neutrino flux. This explains the existence 
of islands. For values of Am^ in between islands, the survival probability 
either increases and decreases for varying Earth-Sun distance, or monotoni- 
cally decreases. The exact location of the islands and their widths can only 
be understood by analyzing the fit procedure, in particular the minimization 
of with respect to the background rate and the incoming neutrino flux. 
Note that there are no "islands" above Am^ ^ 2.5 x 10"^'' eV^. This is be- 
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cause when the oscillation length is small enough (or Am^ large enough) , the 
survival probability cannot only increase for increasing Earth-Sun distance, 
but necessarily reaches a maximum before the aphelion, and then decreases, 
independent of what the survival probability at the perihelion is. This situ- 
ation is qualitatively different from the low point. 

In the case of the high point: initially the survival probability is close to 
unity, decreases sharply as the Earth moves further from the Sun, and then 
grows rapidly, reaching a maximum when the Earth is close to its aphelion, 
because the oscillation length is smaller than AL. In this case, little vari- 
ations in the oscillation length, i.e. Am^, produce big qualitative changes, 
including the position and number of maxima and minima. There is still a 
small ambiguity {i.e. two "islands") in determining Am^ for the high point. 
This happens when the oscillation length is such that the minimum of the 
survival probability happens in March/October and the survival probability 
is large enough at the perihelion and the aphelion. The fact that the absolute 
values of the number of recoil electrons detected are different is taken care 
of by varying the signal and the background. 

In conclusion, if Nature chose neutrino oscillation parameters such that 
sin^2^ is large and Am^ ^ few x 10~^° eV^, Borexino should be able to 
measure these parameters independent of the SSM and any knowledge of 
the number of background events, with good precision (especially in Am^). 
If Am^ ^ several x 10~^^ eV^, the determination of oscillation parameters 
is not as precise. Better precision can be achieved at KamLAND, but the 
ambiguity of solutions in the "low" Am^ region still remains. 

5 Exclusion of Vacuum Oscillations 

In this section, we address the issue of what the experiments can conclude 
about vacuum oscillations if no discrepancy from the normal seasonal varia- 
tion effect is detected. In this case, one may be able to measure the incoming 
neutrino flux, as outlined in Sec. 0. Two distinct possibilities will be consid- 
ered: (1) the measured flux is consistent with the SSM prediction; (2) the 
measured flux is suppressed with respect to the SSM prediction. 

In the first case, one would be inclined to trust the SSM prediction of the 
''Be neutrino flux and use it in the analysis to exclude vacuum oscillations. 
This will be discussed in Sec. |5.1| . On the other hand, in the second case, 
it is not clear if the reduced flux is due to MSW neutrino oscillations, an 
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incorrect SSM prediction of the neutrino flux, etc. This will be discussed in 
Sec. 



5.1 If the Flux is Consistent with the SSM Prediction 

We simulate "data" consistent with the SSM and the expected number of 
background events. The relevant numbers are quoted in Sec. 0. The "data" 



are binned into months (12 bins per year), and are illustrated in Fig. |TT 
assuming three years of Borexino running. We then fit to the "data" an- 
nual distributions that include neutrino oscillations for a given choice of 
(sin^ 2^, Am^), plus a constant background. The background rate and the 
incoming neutrino flux may be allowed to float in the fit, constrained to a 
positive number. 

It is important to note that this is the opposite of what was done in 
Sec. ^, where the sensitivity of Borexino and KamLAND to vacuum oscilla- 
tions was studied. There, the simulated "data" were consistent with vacuum 
oscillations, and one tried to fit a non-oscillation prediction to the "data" by 
varying the incoming flux and/or the background. Here, the "data" are con- 
sistent with no oscillations, and one tries to fit the "data" with a prediction 
which includes the effect of neutrino oscillations for fixed (sin^ 26, Am^), by 
varying the incoming flux and/or the background. If both the background 
and the incoming flux are fixed, i.e. not allowed to vary in the fit procedure, 
the exclusion and the sensitivity regions are the same. On the other hand, 
if both the background rate and the incoming flux are allowed to float, the 
exclusion region is expected to be smaller than the sensitivity region pre- 
sented in Sec. |^, especially in the region Am^ < 10^^° eV^. This is due to 
the fact that a large number of points in the parameter space yield an annual 
variation of the flux which is much larger than 7%, but agrees with the 
shape of the normal seasonal variation. If in the fit procedure the signal is 
scaled down to reduce the amplitude of the variation and the background 
scaled up to increase the number of events, a good fit to the no oscillation 
case can be attained. 

Fig. |l^ shows, for three years of Borexino and KamLAND running, the re- 
gion of the (sin^ 26, Am^) parameter space excluded at 95% CL, if one allows 
the solar neutrino flux and the background rate to float within the positive 
numbers (in white), and if one assumes the solar neutrino flux calculated in 
the SSM within theoretical errors (in light plus white). 

A few comments are in order. First, one notices that the KamLAND 
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Figure 13: Region of the two neutrino oscillation parameter space excluded 
in the case of no neutrino oscillations if one assumes no knowledge of the 
background and no knowledge of the SSM (white) or knowledge of the SSM 
(light+white), after 3 years of Borexino (right) and KamLAND (left) run- 
ning. 



exclusion region is larger than the one excluded by Borexino. This is, of 
course, expected because of KamLAND's larger fiducial volume and therefore 
higher statistics. Second, when the solar neutrino flux is allowed to vary 
in the fit, the excluded region of the parameter space shrinks, as expected 
and discussed earlier. Third, one can safely claim that, if no discrepancies 
are detected in the seasonal variation spectrum, the "large" Am^ (several 
xlO~^° eV^) set of vacuum solutions (see Figs. |^ and |10D will be excluded, 
even at Borexino. Even when no knowledge of the incoming neutrino flux 
is used, a reasonable portion of the "small" Am^ (several xlO~^^ eV^) set 
of solutions is also excluded. When one assumes knowledge of the incoming 
neutrino flux, the entire allowed region is excluded. 

If the background rate is larger than expected, the excluded region dimin- 
ishes accordingly. This is because when the constant background is enhanced 
with respect to the oscillation signal it is easier to achieve a reasonable 
for the fit even when the seasonal variations due to vacuum oscillations are 
significantly different from the no-oscillation case. In particular, when the 
background rate is large enough that the seasonal distribution of the data is 
statistically consistent with a fiat one, a reasonable for the fit can always 
be achieved simply by scaling the signal to zero and scaling up the back- 
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ground appropriately. Explicitly, after three years of Borexino (KamLAND) 
running the exclusion region vanishes if the background rate is ~8 (40) times 
larger than anticipated, when both the background rate and the incoming 
neutrino flux are allowed to float in the fit or ~500 (3000) times larger than 
anticipated when one assumes the neutrino flux predicted by the SSM. 

5.2 If There is an Overall Suppression of the Flux 

If there is an overall, i.e., time-independent suppression of the flux (which is 
the case for the MSW solutions) , the way to proceed towards excluding part 
of the vacuum oscillation parameter space is less clear. This is because such 
an experimental result neither agrees with the SSM prediction nor does it 
represent any "smoking gun" signature for neutrino oscillations, as is the case 
of anomalous seasonal variations. One does not know if the SSM prediction of 
the flux is simply wrong, or if there are neutrino oscillations consistent with 
one of the MSW solutions or both. Anyway, it is clear that (in general) the 
incoming neutrino flux should be considered unknown in the data analysis. 

The most conservative option is to follow the same analysis done in the 
previous subsection, and allow both the incoming neutrino flux and the back- 
ground rate to float in the fit. In this case, the excluded region of the 
two-neutrino oscillation parameter space is reduced significantly, and may 
completely disappear. This is because when the number of signal events is 
reduced the annual distribution is closer to fiat and a good fit is obtained even 
when the would-be annual variations are very different. This is very similar 
to what was previously discussed at the end of the last subsection, where we 
discussed what happens if the background rate turns out to be much larger 
than anticipated. Explicitly, after three years of Borexino running and a 
signal rate which is 21.3% of the SSM prediction (as one would obtain in the 
case of the small angle MSW solution), Borexino is unable to exclude any 
portion of the vacuum oscillation parameter space, while KamLAND can still 
exclude about one half of the "high" and "low" Am^ preferred regions. If 
the background rate can be estimated by other means with 10% uncertainty, 
Borexino and KamLAND will be able to exclude the entire "high" Am^ region 
and a significant portion of the "low" Am^ region. 

In order to go beyond the most conservative analysis discussed above, 
one would have to look at the overall situation of the solar neutrino puzzle 
at the time of the data analysis. It is likely that one will be able to do 
much better. For example, solar neutrino oscillations might have already 
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been established by the SNO experiment and perhaps it is reasonable 
to assume the incoming solar neutrino flux predicted by the SSM. Then it 
would be possible to exclude a region of the parameter space as large as 
the one in Sec. |5.1| where one assumes the SSM flux. Another possibility 
is that Super-Kamiokande or SNO rules out the small angle MSW solution 



by studying the distortions of the electron energy spectrum |3T[|, and a 
large suppression of the ^Be solar neutrino flux would indicate that there is 
something wrong with the SSM. In this case, it is not clear how to proceed. 
We do not go into further discussions on all logical possibilities. 



6 Conclusions 

We have studied possible uses of the seasonal variation of the ''Be solar neu- 
trino flux at Borexino and KamLAND. Our results can be summarized as 
follows. Once the experiments accumulate enough data to see seasonal vari- 
ations, the flrst step will be to determine if the observed pattern is consistent 
with the normal flux suppression. If a discrepancy is found, it will be 
a sign of vacuum oscillations. In this case, the seasonal variation of the data 
can be used to determine the oscillation parameters sin^ 26 and Am^. On the 
other hand, if the data are consistent with the normal pattern, the amplitude 
of the variation can be used to measure the ^Be solar neutrino flux and to 
exclude a signiflcant portion of the vacuum oscillation parameter space. 

If the observed seasonal variations are consistent with the normal 
flux suppression, one can use the amplitude of the variation to determine 
what fraction of the observed recoil electrons are induced by the neutrinos 
coming from the Sun. This method is limited by statistics, and the accuracy 
is worse when the "^Be solar neutrino flux is suppressed, as in the case of 
the small angle MSW solution. In fact, in Sec. ^ we found that in the case 
of a large suppression only KamLAND should be able to perform such a 
measurement, after 3 years of data taking. It is important to emphasize that 
we assumed the oscillation of electron neutrinos into other active flavors. 
In the case of oscillations into sterile neutrinos, the ^Be solar neutrino flux 
might be almost absent, and in this case neither Borexino nor KamLAND 
are able to perform a measurement of the flux using this technique. 

An important advantage of this technique is that it does not require a 
separate estimate of the background rate, which may be a very difficult task. 
If the background rate can be reliably measured by some other means, one 
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can obtain another measurement of the neutrino flux. In this case, the two 
results can then be compared for consistency, thus making the flnal result on 
the ^Be neutrino flux much more trustworthy. 

We also studied in great detail the effect of vacuum neutrino oscillations 
on seasonal variations. Our analysis shows that the outlook for discovering 
vacuum oscillations at both Borexino and KamLAND is very favorable. A 
very important flnding in Sec. |^ is that the experiments may detect a devia- 
tion from the normal pattern of seasonal variations even without relying on 
the SSM prediction of the incoming neutrino flux or estimate of the back- 
ground rate. The analysis would consist of trying to flt the observed data 
with the normal pattern, treating the incoming neutrino flux and the 

background rate as free parameters. With this technique, after three years of 
running Borexino should detect anomalous seasonal variations for almost all 
values of (sin^ 26', Am^) preferred by the analysis of the neutrino flux data 
from Homestake, GALLEX, SAGE, and Super-Kamiokande, as illustrated 
in Fig. p. The sensitivity region should be larger at KamLAND (Fig. 10). 
Results obtained in this way would be very robust. Both experiments are 
sensitive to an even larger portion of the parameter space if the background 
rate can be reliably estimated by auxiliary measurements. 

If anomalous seasonal variations are discovered, the data can be used to 
measure the oscillation parameters (sin^ 20, Am^). This issue was studied in 
Sec. ^. It was found that for Am? ^ eV^ the experiments will be able to 
determine Am? with good precision. At the same time, for Am^ ^ 10~^° eV^ 
there would be many "candidate islands" in the (sin^ 26, Am^) plane, and it 
will not be easy to resolve the ambiguity. 

On the other hand, the absence of anomalous seasonal variations of the 
^Be solar neutrino flux data can be used to exclude regions of the vacuum 
oscillation parameter space. In Sec. ^ we presented the exclusion plots for 
both Borexino and KamLAND, after three years of running. An important 
lesson from that section is that in order to exclude a large portion of the 
preferred region, the experiments will need to either measure the background 
rate or rely on the SSM prediction for the neutrino flux. In the absence of 
both, the results are rather weak. This is to be contrasted with the situation 
in Sec. ^. 

It is important to keep in mind that the simulated "data" is most of the 
time based on the SSM prediction for the ''Be solar neutrino flux and the 
anticipated number of background events at Borexino and KamLAND. Our 
numerical results, therefore, even in the cases when we do not use the knowl- 



31 



edge of the incoming neutrino flux or the background rate at the analysis 
stage, are not to be regarded as SSM and background rate independent. We 
would like to draw attention to our comments at the end of Sees. and |^ on 
how our results might change if these inputs are changed. We also assume 
only statistical errors in the data analysis, neglecting systematic uncertain- 
ties due to the lack of knowledge in the seasonal variation of the background 
rate. The inclusion of such effects is beyond the scope of this paper. 

Overall our results indicate that the future Borexino results can lead to 
significant progress towards solving the solar neutrino puzzle. Furthermore, 
if KamLAND is also able to study solar neutrinos, one would have access to 
a larger data set, and more powerful results can be obtained. 
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A Analysis 

In the analyses in Sec. H, |, and |^, we are interested in the capability of 
an "average" experiment. It is possible to simulate "data" with statistical 
fluctuations included, but then the value of would vary slightly between 
different repetitions of the same simulation. A better approach is to find an 
expression for "averaged" over many simulations. As we show below, aver- 
aging over statistical fluctuations simply leads to the inclusion of a constant 
term in the definition of x^- 

Suppose we have some solar neutrino data binned into A^bins bins. Let 
the average expected value in the ith bin be di with corresponding random 
fluctuation Adi. Suppose we want to fit this data with a function /, which 
can depend on two parameters: the signal s and the background b. Then the 
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of the fit can be defined as follows: 



2, ^-(rf, + Ad,-/,(s,6))2 



(A.l) 



where 0"^. = ^/di + Adi. Because, in the case of interest, the number of events 
per bin is sufficiently large, we can approximately set a^. ~ v^-Q 

First consider the case when s and b are fixed numbers. The average 
value of the one would obtain after simulating the data many times is 



X 



E 



Using (Adi) = 0, {{AdiY) = di, we find 

{d.-f^?' 



{Ad,f ^ 2Ad,{d,-U) ^ {d, - fif 



X 



A^bins 

E 



1 + 



di 



di 



i 



(A.2) 



di 



(A.3) 



Therefore, in this simplest case it is enough to use the average values di and 
the number of bins to compute (x^)- 

Next, consider the case when f{s,b) = b + g{s) and is minimized with 
respect to b. 



dx\s,b) 
db 



i di 

(rf, + Arf,-(7,(g) )^ 
V di 



? d,. 



(A.4) 



Introducing Ai = {di+ Adi— gi{s)) / di and substituting Eq. (|A.4|) in Eq. ( |A.1|) , 
we obtain 



A mm 



A'bins 

E 

i 

A^bin 



'A^bin 



'A^bin 



Ai 



^ Affair 



A2 a / Afbins A 



A 



'Afbini 



^= Ai 



^Afbin 



d, 

-1 



di 



di 



(A.5) 



"One can easily estimate the resulting relative error in to be of 0{l/y^ (di)). 
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Now plugging back in the definition of Ai, we perform the averaging using 
{Adi) = 0, ((Arf,)2) = d„ and ((At/,) (At/,)) = for z ^ j: 




Nuns + E (A.6) 
-A J-M^V (A.T, 



Substituting Eq. ( [A.6| ) in Eq. ( |A.5|) , we find 



The last two terms are exactly what one would find after minimizing 

b — gi{s)Y /di with respect to 6, and hence in this case random fluctuations 

can be accounted for by replacing iVbins in Eq. ( [A.3| ) by A^bins — 1- 

One can easily show that, if /(s, h)i = b + s-hi and one minimizes with 
respect to s, the effect of random fluctuations is also to substitute Abins — 1 
for Abins in Eq. (|A.3|) . The proof is completely analogous to the case we 
just studied. Moreover, it is straightforward to combine the two results and 
consider minimization with respect to both b and s, in which case one should 
replace Abins in Eq. (|OD by A'bins - 2. 

In general, one should use the number of degrees of freedom A^d.o.f. when 
computing (x^): 

X') = iVd.cf. + E ^^M^. (A.9) 

i 



di 



B Analytic Estimate of the Sensitivity Cutoff 

In Sec. ^ we showed that the sensitivity region for anomalous seasonal vari- 
ations is limited by the finite linewidth of the ''Be line. In this appendix we 
show how one can analytically estimate the location and the shape of the 
sensitivity cutoff. 

As was mentioned in Sec. ^ the true shape of the ''Be line is rather 
complicated, with a Gaussian profile on the low end and an exponential tail 
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on the high end. For the purpose of this estimate we choose to approximate 
the Gaussian part by a sharp cutoff: 



fiE) 



if E < El 



(B.l) 



To determine the fraction of neutrinos reaching the Earth we integrate 
the oscillation probability P{E, L) given by Eq. (|3.1|) over the line profile 
Eq. ( [B.l| ) and divide by the normalization constant N. 



P{L) = - I dEP{E,L)f{E) 
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N 
+ 



sin2 20' 



1 - 



dEe 



-aE+b 
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dE cos 2 



El 



1.27Am2L 



{Eo -E)]e 
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sm^20\ 1 sm^20 



1 



N 



( 1.27 Am^L^ 
xcos (2 {Eo 



El) — arctan 2- 



a2 + (1.27Am2L/Eo2)2 
1.27Am2L\^ 



aEi 



(B.2) 



Since the width of the line is only several keV while Eq = 0.862 MeV, we can 
set Eq — El ~ Eq in the argument of the cosine. Substituting the value of the 
normalization constant = J^^dEe~°'^~^'' = (1/0)6""^^"*"* and introducing 
(j) = arctan(2 x 1 .27 Am'^ L / (aEg)) , we obtain 



P(L) ~ 1 



sin^ 29 



cos 



(2i 



27Am?L 
Eo 



1 + {1.27Am^L/{Eia)y 



(B.3) 



From this equation we can read off the shape of the cutoff. Viewed as 
a function of Am^, for small values of the mixing angle the cutoff profile is 
approximately given by 



sin^ 2^^eutofr oc a/I + {1.27Am^L/{Eia)y. 



(B.4) 



Using the numerical value of a = 0.75 keV ^, obtained by fitting the line pro- 
file in ||2^, we find that sin^ 26'cutoff(^"^^) should increase by ^/2 with respect 
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to the smallest value of sin^ 26'cutofr when Am^ ~ 2.9 x 10~^ eV^. The actual 
number from curve 4 in Fig. |^ is Am^ ~ 1.5 x 10^^ eV^. The actual value 
is smaller, which is expected, because, for the purpose of this estimate, we 
neglected the contribution of the Gaussian part of the line profile, effectively 
making the line narrower. 

One can also estimate the location of the cutoff if the line profile were 
purely Gaussian (curve 3 in Fig. ^. The steps are completely analogous: the 
new normalization constant is A^' = dEe~^^~^°^ I" = y/ncr, and P{L) 
is given by 

P(L)^—J_JEP(E,L)e-'''''' 

N' [\ 2 J J-oo 

+ J_^dE cos \2 {Eo - E)\ e"^ 

= 1-^(1- e-(-^^-^ cos (2^^^^^^) ) . (B.5) 

Thus, the cutoff for this model sets in faster and the profile for small values 
of sin^ 29 is Gaussian. Numerically, sin^ 26'cutofr(Am^) is expected to increase 
by -\/2 with respect to the smallest value of sin^ 2^cutoff when Am^ ~ 4.2 x 
10^^ eV^, which agrees with curve 3 in Fig. 0. 
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